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It is well known that the strong subadditivity theorem is hold for classical system, but it is very 
difficult to prove that it is hold for quantum system. The first proof of this theorem is due to Lieb 
by using the Lieb's theorem. Here we use the conditions obtained in our previous work of matrix 
analysis method to give a new proof of this famous theorem. This new proof is very elementary, 
it only needs to carefully analyse the minimal value of a function. This proof also shows that the 
conditions obtained in our previous work are stronger than the strong subadditivity theorem. 
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I. INTRODUCTION 



Entropy is an important concept not only for physics but also for information science. From the definition of 
entropy, we can get some simple properties of it, such as concavity, continuity property, additivity and subadditivity 
[1]. But some other properties is not so obvious, such as the strong subadditivity theorem (SSA). Among all of the 
properties of entropy, the most famous one is the SSA, and it is very difficult to prove this theorem for quantum 
\ system. The content of this theorem can be expressed as the following: two overlapping subsystem AB and BC, the 
I/"") ■ entropy of their union (ABC) plus the entropy of their intersection (B) does not exceed the sum of the entropies of 
the subsystems (AB and BC) [2], that is 



O. S(p AB c) + S( PB ) < S(pab) + S(pbc)- (1) 

"xT ' 

where S(p) = —Tr(plnp). It is well known that this theorem is true for classical information theory, but to prove this 
' theorem is true for quantum system is very difficult. This theorem is first conjectured to be true for the quantum 
Oh, system by Lanford and Robibson [3]. The first proof of this conjecture is given by Lieb et al. several years later. This 
proof is based on the concave of the function S(p\2) — S(p\) in pi2 [4]. Another proof based on the same fact was 
proposed by Uhlmann [5,6]. 

Recently, quantum information theory attracts more and more attentions for its misterious properties and its 
£7^" potential applications in science and technology [7,8]. The SSA plays an important role in this new field [9] too. 

The fundamental scource in quantum information is entanglement between many particles which can be viewed as 
i— i , the relations between the partial particles. So distinguish whether a set of the partial particles come from a single 
state (Y-representability problem) [10,11] and further to obtain its entanglement property are important in quantum 
information while SSA gives a strong constraint on the partial particles and the whole system. The convenience of 
the SSA is that it has explicit physical and manipulating meaning. So it is a convenient necessary criterion for the 
Y-representability problem. Recently, we use the matrix analysis method [12] to get some necessary conditions for 
the Y-prescntability problem. We find that using these conditions we can get a new proof for the SSA. Our new proof 
is elementary, we need only to use the Lagrange multiplier method and carefully analyse the minimum of a function. 



II. THE THEOREM AND THE PROOF 



There is a density matrix pabc, where the particle A, B and C are in L— dimension, M— dimension and 
N— dimension Hilbert space, respectively. Let {A^g, X^ B , • • • , A^ 1 '}, {A^, A^l, • • • , A^ 1 ^'}, {A^ , A^ , • • • , A^} 
and {XabC , ^abC > ' ' ' > ^abc*^ } are ^ ne eigenvalues of the density matrix pab , P bc > P b and pabc > respectively (where 
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Pab, Pbc and Pb are gotten by tracing the other particles from pabc)i and they are arranged in increasing order. We 
defined vectors X AB = {A^, X%, • • • , A^°}, X BC = {A^, X% ] c , • • • , A^'}, X B = {a£\ A^ 2) , • • • , A^ M) }, X A bc = 

r>(l) x (2) .(LMN)-, , x s _ r^L .(») »(») r IM x (i) 1 X B \(») 

I ABC ABC '"' ABC I duu AB ~ tZ^i— 1 AB 1 2-^i—L + l AB> "' Z^j=L(M-l)+l A AB J ' A BC — 1^=1 A BC 
^2JV x (i) tt^MN x (i) , X AB _ \0) V- 2JV \C?) s-^LMN ,(j) -. , 

Z^i=N+l A BC> "■' Z^i=(M-l)A r +l A AB/' ABC — 1^=1 A ABC 2^j—N+l ABC> Zjj-(LM-1)jV+1 ABC J> dnu 

A abc ={Ej L =i a abc E^l+iAabc E™mjv-i)l + i A abc }• Usm S thc matrix analysis method, we get the 
following two lemmas on the eigenvalues [12]. 

Lemma 1. Using thc notes defined before, we can get the relations between the eigenvalues of pbc, Pab, Pb and 
Pabc as 

Aabc >~ Xab (2.1) 

Aabc >~ Xbc (2.2) 

Afc y \ B (2.3.1) 
Xab >■ Ab (2.3.2) 

Lemma 2. Suppose rank(pABc) — LMN—Ls, rank(psc) — MN—s, rank(pAB) = LM—r and rank(pB) = M—t, 
if r and s satisfy the condition Nr < Ls, there will be 

t<[^] + l, (3) 

where [x] is the maximum integer which is smaller than x. 

The notation y >~ x mean that the vector x is majorized by the vector y. The majorization is defined as the 
following. Let x — {x\, X2, ••• , x n } and y = {y\, y 2 , • • • , y n } are n-dimensional vectors and the elements are arranged 
in increasing order. Then the vector x is majorized by vector y [13], denoted by y >~ x, if for each k (k = 1, 2, • • • n) 
the following inequality is hold 

k k 
i=l i=l 

and the equality is hold when k = n. Under these two Lemmas, we can find that the SSA is hold in the following. 

Theorem There are four normalized vectors X ABC = {X ABC , X ABC , • • • , X ABC N }, X AB = {X AB , X AB , ■ ■ • , X AB }, 
X BC — {Af c , X BC , • • • , Xff N } and X B = {Xf , X B , ■ ■ ■ , X M }, the elements of these vectors are non- negative and ar- 
ranged in increasing and define the vectors X AB P, Xg BC and X^ c , X^ B , which are similar to the vectors in lemma 1. 
If the elements of these vectors satisfy the following conditions 

1. X ABC y X AB ; 

2. X ABC y X BC ; 

3. X BC y X B and X AB >- X B : 

4. Suppose the vector A = X ABC has only Ls zero elements and X BC has s zero elements, and if the vector X AB has 
r zero elements, there are at least [^-j^] + 1 elements of the vector X B are zeroes. If exchange the role of the vector 
X AB and X BC , the similar result must be hold also. 

Thus the following inequality is hold 

S(X) + S(X B ) <S(X BC ) + S(X AB ), (4) 

where 5(A) = J2f^ N (-Xi In A,). 

The proof of the theorem is technical. We use the Lagrange multiplier method to get the minimal value of a function 
under the conditions 1, 2 and 3. Because there are many possible extreme points, we need to find out the minimal 
one. We use some facts to find that when the function gets the minimal value, all of the nonzero elements are equal 
to each other. Then use thc condition 4 to get the minimal value of thc function. 

Proof. At first, we define a function 

LM MN M LMN 

F = J2(-Xf B lnX AB ) + 5>Af c lnAf c ) + £(Af InAf ) + £ (A, In Aj), (5) 

2—1 i—1 i—l i—1 
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If we can prove that the minimum of this function is not less than under the conditions 1-4, the theorem is true. 
So the proof becomes to find the minimal value of a function under some conditions. Obviously, the minimal value of 
this function exists and is finite. Now we use the Lagrange multiplier method to deal with the conditions 1, 2, 3 and 
define a new function 

LM MN M LMN 

2 — 1 l—\ 1—1 l — l 

LM-1 i Ni MN—1 k Lk 

+ E c4(E a / s -E a -*h)+ E #(E A f°-E A i-&) 

i=l j=l j=l k=l j=l 3=1 

M-l Li i M-l Nk k 

+ E °?(E A f - E A / B - 4) + E ^(E A f - E A ? c - 

LMN LM 

+ E 7i(Vn — Aj — zf) + ]T ^(Af + B i - 5 - r?) 

j=0 i=l 
MN M 



+ y: \f c - o + E - A ? - *?) 

i=l i=l 

LM MjV M LMN 

+ ai (J2 Xf B 1) + a 2 (£ Af c - 1) + a 3 (E A f " !) + <*( E A * ~ ^ 



i=l i=l i=l i=l 

where the parameters ct\, (5 2 k , Ui, Vi and Wi, "fj, di are Lagrange multipliers, y\ k , y\ k , zf, rf, sf and if are 
introduced to make the inequalities to be equations. We have used the conditions that X AB ,X B , X BC and A 2 - are 
arranged in increasing order and let Ao = 0. 

Then when G get the minimal value, there must be some constraints on the parameters and variables. First, we 
can get afxji = (i = 1, 2, • ■ ■ , LM — 1; j = 1, 2) and the similar relations between (31 and yjk, 7i and Zi, Ui and ti, Vi 
and Si, Wi and ti. The second, we get the relations between the elements of the vector X AB and A, 

i Ni 

J2 x t B -J2^ = x2 ^ i = 1 ^---' LM - 1 w 

and the similar relations between A and X BC , X AB and A B , X BC and X B . Then the relations between the vector X AB 
can be gotten 

X A + B -X AB =rf,i = l,---,LM-l; (8) 

LM 



£^-1 = 0, 



i=\ 



and the similar relations between the vectors of X BC , X B and A. 

The most important constraints are the equations between the vectors Xf B , Af c , Af and A; 



LM-l M-l 



-\n\ AB - 1 + a i - E (rf+Ui- 1 -Ui + a 1 =0(i = l,2,---,LM). (9.1) 



MN-l M-l 

lnAf c -l+ P}- E Pj +Vi-i-Vi + a 2 = (i = 1,2, MN). (9.2) 

3=1 j=[Vl+ 1 



M-l M-l 



InAf + 1+ E a * + E +w t -i-w l +a 3 = (* = 1, 2, • • • , M). (9.3) 
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LM-\ MN-1 

InAj + 1- a )~ 2 / 3 ]+7i-i-7i + a 4 = (i - 1,2, • • • ,LMJV). (9.4) 

j=[V-]+i M^l+i 

Since the number of the possible cases are so large, it is very difficult to get the solutions directly. We point out 
some useful facts to reduce the possible solutions and to find the minimal value of the function G. 

Fact 1. When the function G get the minimum, suppose that parameters a} and ctg are the nearest nonzero 
parameter act on the elements of vector X AB , if the parameters Ui and ULj are zeroes, all the parameter Uk (i < k < Lj) 
are equal to zeroes. This fact is true for the other parameters Vi, u>i, 7$. 

Proof. Without loss of generality, we only consider the parameter U;. Suppose the fact is not true, there are some 
parameters u p (i < m < p < n < Lj) are not zeroes. For simplicity, we suppose there are no more nonzero parameters 
aj and o?-. Then we get the conditions from (9.1) 

- In X AB - 1 + a] - a) + - u m + <n = 0, 
- In A^+! - 1 + a} - Oj + u rn - u rn+1 + ai = 0, 

(10) 

- In \ AB - 1 + a] - a? + u n -i - u n + ai = 0, 

- InA^fj - 1 + a] - a 2 j + u n - + a x = 0, 

where we have used the conditions that the parameters u m -i and u n+ i are zeros. Since the parameters u p (i < m < 
p < n < Lj) are nonzero, then we get = X^+i = ■ ■ ■ = \ AB = ^n+i- So we have the relations between these 
nonzero parameters 

-u m =u m - u m+ i = ■■■ = u n -i -u n = u n , (11) 

that is, u n = (n — m + l)u m = —u m . So u m = 0, then all of the parameters u p (i<m<p<n< Lj) arc zeroes, 
which is inconsistent with our suppose. So this fact is true. QED. 

Since the fact 1, the parameters Uk affect the result only when there are some nonzero parameter a} or ctg make 
k = i or k = Lj . For this situation, we have the following fact. 

Fact 2. When the function G get the minimum, if there are a set of parameters Uk (m < k < n) are nonzero and 
there are some nonzero parameters a\ make m < i < n, This situation is equal to the situation where the parameters 
Uk (m < k < n) and a\ are zeroes, but two new parameters aj n _ 1 and a\ +1 should be added, and the parameters 7, 
(0 < i < LMN) should be adjusted. 

Proof. Without loss of generality, we suppose only the nonzero parameter a] satisfy the condition m < i < n. For 
simplicity, we suppose there is no other nonzero parameters act on the eigenvalues of Pab- Since the parameter a] 
are nonzero, then 

i Ni 

£Af*-5>=0; (12) 
and the parameters Uk (m < k < n) are nonzero, we get 

\AB _ \AB _ _ \AB _ \ AB /i o\ 

A m — A rn+1 — " ' — A n — A n+1- 

Since we have the condition J2]=\ A f B — J2f=i~ 1 ' > A j> together with the equation (12), we get \f B < 

T,f=N(i-i)+i x j- 0n thc othcr hand ' x t+i > ESVj. that is > EjL'jv(i-i)+i X o > ES+i V Because of the 
condition A^^.^+i < Ajv(i-i)+2 < • ■ • < A N(i+i), wc get the equation Ajv(i-i)+i = Ajv(i-i)+i = ■ • ■ = Ajv^+i). So we 

get \? B = Ef=V(i-i) + i A j and A^ = A j, that is 

i-l JV(i-l) i+i JV(i+l) 

5>f - E A^O^A-- 2 A,=0. (14) 

i=i j=i i=i i=i 

Continue to use this method we can get 
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m-\ N(m-l) n +l N(n+1) 

E X t B - E A,=0,£Af - £ Ai = 0. (15) 

i=i i=i i=i i=i 

From these equations, we can find this is just as there are two nonzero parameters a} n _ 1 and a* +1 , and the nonzero 
parameter a\ have no effect in this case. From the constraints on \£ B , there wll be 

-lnA^ B - 1 + a\ +ai = (1 < k < m - 1) , (16.1) 

- In \£ B - 1 + a] + u k -i -u k + ai = (m<k<i), (16.2) 

- In \£ B - 1 + u fe _i - u fe + ai = (i + l<fe<n + l), (16.3) 

- In \£ B - 1 + ai = (n + 2 < k < LM) . (16.4) 

Since the elements X AB (m < k < n + 1) are equal to each other, then u m — u m — u m+1 = ■ ■ ■ = — iti = a, 
Ui — u i+ i = u i+ i — u i+2 = ■ ■ ■ = u n = (3 and (3 = a + aj. If we let the nonzero parameters a] n _ l = —a and a\ +1 — [3, 
the equations are the same. Now we consider the effect of this substitution on the vector A. For simplicity, we suppose 
also that there are only the nonzero parameter a\ act on the vector A. Then the equations are 

\n\ k + l-a\+ lk _ l - lk + a i = (1 < k < Ni) , (17.1) 
In A fe + 1 + 7fe _! - 7fc + a 4 = (Ni + 1 < k < LMN) . (17.2) 

Insert the parameters a} n _ 1 and a^ +1 into the equations, we can find that 

In A fe + 1 - - al +l + 7 fe_! - 7fe + a 4 = (1 < k < Ni) , (18.1) 

mA fe + l-a£ +1 +7fe_ 1 -7fe + a 4 = (Ni + 1 < k < N(n + 1)) (18.2) 
In Afe + 1 + 7fe_! -7I + a 4 = (N(n + 1) + 1 < k < LMN) (18.3) 

where 7 fe (Ni + 1 < k < LMN) are the new parameters to make the equations are the same as the equation (17). 
This is just as the situation that the parameters aj re _ 1 and a^+i are nonzero, and the parameter 7fe is adjusted. QED 
This fact is also true for the parameters /3j. This fact tell us that any solution found in the former situation can 
be found in the later case. In the following, we always suppose we have already done this change. After making 
these changes there is no nonzero parameters a] or /3j make the parameter Ui(vj) nonzero and the parameters ji are 
substituted by 7-. 

Fact 3. When the function G get the minium, if i and j are the nearest indexes to make the equations 
J2k=i \^ B ( BC ^ — ^fe^i^ and J2k=i \ BC ^ AB ^ — J2k=i ^ A fe to b e hold, the elements of the vector A between 
Ni and Lj are equal to each other. 

Proof. We suppose this conclusion is not true, without loss of generality, let Lj > Ni. Then there are some elements 
satisfy the following conditions \ L j = ^Lj-i = • • • = \ = A fc > A a = X Ni = A^i+i = ■ ■ ■ = X q , for simplification, we 
suppose that Lj — p > q — Ni. If we define the following parameters Aj and A m as Y^k=i+i ~ Y^k=Ni+i ^fe = A; 
(i + 1 < I < , [x] is the maximal integer which is smaller than x) and J2k=Nm+i ^fe ~ J2l=m ° = 

{[^jf\ w can find that all of these parameters are more than zero. Then we take out the minimal number 

from and , we denote it by A, obviously it is more than zero. Now we change the element A2 by A2 — ^f_ p 
and Ai by A 4 + - A Ni where the parameter A = (q — Ni)A. After these substitution, the new elements of the vector 

A satisfy all of the conditions. The entropy of the vector A is larger than the entropy of the vector A and the entropy 
of the other vector is invariable. So the function G for the new vector is smaller than the former which is inconsistent 
with the suppose. QED. 

This fact is also true for the vector X B . Since we have the fact 3, then we want to know how many nonzero 
parameters a] and /3j in the section where all of the elements are the same. We have the following fact 

Fact 4. When the function G get the minimum, there is no nonzero parameters a} and (3j in the section where all 
of the elements of vector A are the same except for the edge parameters. 

Proof. We first point out that there are at most four nonzero parameters a] or (3j in the section where all of the 
elements of vector A are equal to each other if the conclusion is not true. If this assert is not true, there are at least 
five nonzero parameters act on the section where all of the elements of the vector A are the same. So at least three 
of them (such as a] (I = i,j,k---) or f3j (I = i,j,k---)) are act on the same vector. Without loss of generality, 
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we suppose there are three nonzero parameters a* (I — i,j,k). Since the elements Ajvi+i — = ■ ■ ■ = Ajvj = 

•^JVj+i = Ajvj+2 = • • • = X^k = A a - We have A^f > NX a and X AB < NX a , since A^ < X AB , then all of the elements 
X AB (i + 1 < I < k) are equal to each other. Because we have already done the changes in the fact 2, and use the fact 
1, we find all of the parameters Uk (i + 1 < I < k) are zeroes. Further more, the parameters are zeroes too. So the 
number of the nonzero parameters is no more than four, and they divide the section where all the elements are equal 
into three smaller sections. 

Now we only need to prove the case that less than five parameters are also zeroes. If these parameters are nonzero 
and set on the vectors as figure 1, which makes the function G get the minimum, wc take some sufficient small value 
A from the elements of the first section to the third section. At the same time, A must be taken from the left side 
section of the parameters k and j to the right side section. Using the same method of the proof of the fact 3, if the 
A is sufficient small, all the conditions will be satisfied. From the following calculating, we can find that through this 
manipulation the function G is smaller which is inconsistent with the minimal suppose. 

Let the elements of the vectors before the manipulating are Aa^+i = Aa^ + 2 = • •• = X^k = A a , \Nj+i = = 

— \ — \ ■ \AB _ \AB _ _ \AB — \AB \AB _ \AB _ _ \AB — \AB. \BC _ \BC _ _ \BC - 
■■■ — ALl = Ab, Aj — Aj-i — • ■ ■ — A s — A a , A j+1 — A j+2 — ■ ■ ■ — A t = A fe , A k — A k _ 1 — ■ • • — A u = 

A a C > A k+i — A k+2 = ■ ■ ■ = X BC = Af c . After the manipulate, the new elements are A a = A a — Lk A Ni , X b = A a + u A N . ; 
K B ' = ^ B K B ' = K B + A a BC ' = Af c Af c ' = Af c ^ and the other elements are the 

same as before. Since A is sufficient small, we can expand the function ln(A+ ^) = In A+ in the first order. Using 
this formula, we can calculate the difference of the function G between these two vectors. 

G -G= -{Lk - Ni)\ a \n\ a - {LI - Nj)\ b \n\ b + {Lk - Ni)\ a \n\ a + {LI - Nj)\ b \n\ b 
+{j s + l)\ AB In X AB + {t- j)X AB In X AB {j s + 1)X AB ' In X AB ' 
-{t- j)X AB ' In X AB ' + {k-u + 1) Af c In Af c + {v - k)X BC In Af c 
-{k-u+ 1)X BC ' In Af c ' -{v- k)X BC ' In X BC ' 

\^\AB\BC 

= Aln b ° a (19) 

^ m \ \AB\BC [li) > 
A a A b A b 

Since A Q = Ab and X AB < X AB , Af c < Af c , there will be G' — G < 0. This is inconsistent with the suppose that 
the function G get the minimum. QED 

Fact 5: When the function G get the minimum, there are at most one a\ and one a* are nonzero and the elements 
A fc = {k < Ni or k < Lj), X AB = {k < i) , Af c = {k < j) . 

The proof of this fact is similar to the proof of the second part of the fact 4. If there is another nonzero parameter, 
we can take some small value from the left of this parameter to the right of it to make the value of the function G 
smaller, which is inconsistent with the minimal suppose of the function G. This fact means that all of the nonzero 
elements of the vector A are equal to each other. If there is no nonzero parameter act on the vector A, that is, all 
of the parameters a] and (3j are zeroes, then all of the elements of the vector A are L ^ N , all of the elements of the 
vector X AB are j^j , all of the elements of the vector X BC are jjjj , all of the elements of the vector X B are jj . Now 
the value of the function G is zero. If there is only one parameter (such as a}) is nonzero, we have the following fact. 

Fact 6: When the function G get the minimum and there is only one parameter a] (/3j) is nonzero, then all of the 
nonzero elements of the vector X BC , X AB and A s are equal to each other. 

Proof: Without loss of generality, we suppose the nonzero parameter is a] . The nonzero elements of the vector 
X AB is equal to each other. We can get this result by only using the inequality between the elements of the vector X AB 
and A. We focus on the other part of the fact. Since the nonzero elements of the vector X AB are the same, all of the 
parameters af are zero. Now we only consider the parameters Suppose the nonzero parameters @P ,(3f 2 ,---, (3f k 
are set as the figure II. From the constraints of the elements of the vector X B and X BC in equations (9) 

M-l 

-lnAf c -I- Pj+ v i-i -Vi + a 2 = (i = 1, 2, • • • , MTV), (20.1) 

M-l 

InAf + 1+ & ~Wi+a 3 = (i = 1,2, • ■ ■ ,M). (20.2) 

3=i 

Then we find the elements of these vectors can be divided into several groups, in each group the elements are equal 
to each other, that is, 
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\BC _ \BC _ _ \BC — /-BC 

A Nh+ Pl +l — A Ni!+ Pl +2 — • • • — A MN — 1, , 

\BC _ \BC _ _ \BC _ /-BC 

A Nn- qi — A Ni!- qi + l — • • • — A A r il+pi = (, 01 

,BC _ \BC _ _ \BC — /-BC 



\B<J _ \B(J _ _ lift,' _ ABU (m\ 

A Afi 2 +p 2 + l ~~ A Wi 2 +P2+2 — ' " — A Nn- qi -l — Si ! l zl J 



and 



— A 2 — • • • — A Ni k = Sfc 1 



A ii + 1 — A u+2 — — A M — SO 

\ B - \ B - — \ B 
\ 2 + l — A i 2 +2 — — A n 



A f 2 +i — A F 2 +2 — ' ' ' — A n — Ci i ( 22 ) 



A l — A 2 — ' ' ' — A f k — Ck ■ 

We must note that all of the parameters Wi which act on the vector A B are zeroes. At first, if all of the indexes ij 
satisfy Wi j — 0, using the fact 1, all of the parameters are zero. The second, if there are some indexes (such as ij) 
make the parameter Wi- to be nonzero. Because the elements in the same section are equal to each other for the fact 
3, we get Cf_i = Cf • Because of Cf_i > N( BC l3 and (f < N( BC l3 , then ( B = .Y s f So if we let I = [f ] and 
m= we can get the inequality CjLi > (Pj -mXfSij + (N- Pj + m)CjL C i and Cf < (Vj-l)Cf-i,j + (N -qj+l)Cf C ■ 
Since Cf-i ^ C^i j > C ' we can S e ^ that C/^i = Cf^i j — Cf C ■ Now we can get the conclusion by using the fact 1, 
that all of the parameters Vk (ij — q% } < k < ij — Pij) and pf, are zeros. So the second situation can be reduced to 
the first situation. So The constraints of the vectors X B and X AB are reduced to 

M-l 

-lnAf c -l- -Vi + a 2 = (i = 1, 2, • • • , MTV), (23.1) 

j=[V-]+i 

M-l 

In Af + 1 + ^ /?? + a 3 = (i = 1, 2, • • • , M). (23.2) 

If let -^r = x^-7mf = Xi-i.i and 7%- = % (* = 0,2, •••,&), we can get Xi = % and Xi-i,i = Xi-iXT where 
Wi + &i = 1 and = p .+ g . , = ■ From these definition, we find that all of the parameters \i an d Xi-i,i 
are in the section [0, 1]. Since we have the conditions YlfJlm+i A f° = Y^\=i+i A f j tncn we can § et tne equations 

Z ^' = JV 'j+ 1 ' _ Z ^' = 'j + 1 ' That :„ 

yrmv - yrn ^g-- mat is 

[N(ij - i j+l )-p j+1 - qj} X j +Pj+iX°i+iX? +1 + <ljX°i 3 X^-i _ (ij - ij+i)Xj 



MN — Nil - Pi + PiXi 1 M-ii 

So we can get the equations 

(ij i j+ i)pi(l - xD - (M - h)\p j+ i(l (*^r^) + gj (l - (^piDl- (24) 

Xj Xj 

If there is a parameter WMi m+1 = 0, then the mth equations in equations (24) has no item which is including p m +i- 
According to the number of the parameters which make WNi l+1 — (1 < I < k), we can divide the elements of these 
vectors into some sections, the last equation of this section has no item which contains p. We can only point out that 
the parameter Ujvi„ must be zero where the parameter n satisfy the condition rj n = and r/ n _i > 0. Or the condition 
4 will not be satisfied. So we always can sum up all of the equations in the same section to get 

(h i a+ i)pi(i - x ?) = (M - n)[qi(i (— D + j^( Pl + qi)(i ^(X^r <n(—D] (25) 

Xi Xi-i Xi 
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where the parameter s means that the parameter Vi s+1 = 0. We first focus on the lhs. of the equation (25), and 
obviously, it is non-negative. Then we consider the rhs. of this equation, there is a function f(x) = xa 1 ^ x + (1 — x)a~ x . 
The value of this function is not more than 1. Then the rhs. is non-positive. To make the equation to be hold, the 
two sides of the equation must be zero. That is \i = 1 an d = or a = 1. For each section, we can get the same 

conditions which imply that all of the nonzero elements of the vector X B and X BC are equal to each other. QED 

For the case there are two nonzero parameters a\ and /3j , using the similar method before and notice the condition 
4, we can get the same result that all of the nonzero elements are equal to each other. 

For the facts proved before, we can get the conclusion that when the function G get the minimum, all of the nonzero 
elements of the vectors X AB , X BC , \ B and A are equal to each other. Using the condition 4, we can calculate that the 
minimum of the function G is not less than zero. This is the end of the proof of the theorem. QED 

Since the Lemma 1 and Lemma 2, the theorem imply that the SSA is hold. This method can be used to prove 
some other entropy properties between the partial density matrix and the multipartite density matrix, Such as the 
inequality S(pab) < S(p A ) + S(p B )- 



III. CONCLUSION 



In this paper we give a new elementary proof of the SSA which is an important property of the entropy for classical 
information and quantum information. The proof is dependent on the analysis of the minimal value of a function 
under some conditions. This proof also show that the conditions in our previous work [12] are stronger than the SSA. 
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Figure caption. 

Figure 1. In this figure, the parameters a} , a j , ft\ , ft] and , o? q , , /3p are nonzero. The eigenvalues \ v (i < v < j) are 
equal to each other. Theses nonzero parameters divide the eigenvalues between i and j into three sections. We take sufficient 
small value A from the first section to the third section. And the same time, we must take the same value from the left 
section of the parameter k and j to the right section in the eigenvalue X BC &nd X BC , respectively. The bold line means 
that the eigenvalues in the line are the same. 
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Figure 2. In this figure, the parameters ffi , /3^ 2 , • • ■ , 0f are zero. The bold line means that all of the eigenvalues lie in the 
line are equal to each other for the nonzero parameters Vi. 
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